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Abstract 


The  eastern  United  States  continental  margin  profile  Is  relatively 
uniform  throughout  the  Middle  Atlantic  Bight,  but  in  the  South  Atlantic 
Bight  it  bifurcates  Into  an  Inner  and  outer  slope  region.  Coastal  tide 
gage  records  indicate  that  sea  level  oscillations  with  periods  longer  than 
one  week  can  propagate  southward  as  continental  shelf  waves  in  both  Bights 
thereby  providing  a coupling  mechanism  between  the  Bights.  However, 
several-day  period  motions  appear  to  be  confined  to  the  South  Atlantic 
Bight  and  may  result  from  backscattering  of  long  wave  energy  by  the  vari- 
able topography  and  the  Gulf  Stream.  The  coastal  sea  level  phase  data  for 
the  several -day  period  motions  is  not  easily  attributable  to  a monochro- 
matic propagating  wave;  rather,  it  appears  that  wave  group  properties  may 
lead  to  a more  consistent  explanation  of  the  phases.  Cross-shelf  and 
longshelf  wind  stress  components  were  both  strongly  coupled  to  sea  level 
fluctuations  for  long  periods;  short  period  motions  were  more  closely  asso- 
ciated with  dynamic  responses  to  atmospheric  pressure  fluctuations. 


1.  Introduction 


The  eastern  United  States  coastal  and  continental  shelf  zone  may  be 
conveniently  divided  Into  two  major  embayments  or  bights  (Fig.  1).  The 
South  Atlantic  Bight  l SAB)  extends  from  the  Florida  Keys  to  Cape  Hatteras, 
and  the  Middle  Atlantic  Bight  (MAB)  extends  from  Cape  Hatteras  to  at  least 
as  far  north  as  the  Gulf  of  Maine.  Cape  Hatteras  provides  a salient 
coastal  demarcation  between  the  Bights. 

The  "boundary"  between  the  SAB  and  MAB  Is  much  less  distinct  when  one 
considers  the  shelf  and  slope  topography.  The  200  m and  2000  m Isobaths 
in  Fig.  1,  which  approximately  bound  the  most  steeply  sloping  shelf  zone 
in  the  MAB,  have  considerably  less  curvature  than  the  Cape  Hatteras  coast- 
line in  the  Bight  transition  region.  Since  the  continental  shelf  and  slope 
topography  can  be  expected  to  play  an  Important  role  In  determining  low 
frequency  water  movements  through  angular  momentum  (vorticity)  constraints, 
the  Bights  may  in  fact  be  quite  strongly  coupled. 

The  Bights  differ  in  several  dynamically  significant  ways.  As  sug- 
gested by  the  isobaths  in  Fig.  1,  the  continental  slope  Is  a single  coher- 
ent feature  throughout  the  MAB  (and  possibly  farther  north  to  the  Grand 
Banks),  but  In  the  SAB  it  bifurcates  into  an  upper  slope  that  follows  the 
coast  and  a lower  slope  that  turns  offshore  and  becomes  the  Blake  Escarp- 
ment. The  bottom  profiles  along  lines  roughly  perpendicular  to  the  200  m 
isobath  are  shown  In  Fig.  2 for  the  four  Carolina  Capes  and  for  the  Chesa- 
peake and  Delaware  Bay  areas.  The  bottom  slope  changes  by  a factor  of  2 or 
3 between  Capes  Hatteras  and  Romain,  but  relatively  little  over  the  entire 
extent  of  the  MAB.  The  inshore  Gulf  Stream  front  generally  follows  the 
upper  SAB  slope,  and  the  n turns  seaward  at  Cape  Hatteras.  The  velocity 
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Map  of  the  eastern  United  States  continental  margin,  showing  the 
Middle  and  South  Atlantic  Bight  sea  level  and  meteorological 
observation  station  locations.  The  steepest  continental  slope 
zones  are  shown  by  the  200  m and  2000  m isobaths. 
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Topographic  profiles  of  the  continental  shelf.  Profiles  are 
shown  for  the  four  Carolina  Capes  in  the  SAB  and  for  the 
Chesapeake-Del aware  Bay  region,  which  is  typical  of  the  entire 
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structure  of  the  Stream  perpendicular  to  a section  off  Cape  Fear  averaged 
over  two  summer  months  is  shown  in  Fig.  3 (from  Richardson,  Schmitz  and 
Niiler,  1969).  The  maximum  horizontal  shear  of  the  mean  Stream  approaches 
the  magnitude  of  the  local  vertical  component  of  the  earth’s  angular  rota- 
tion rate  (f);  thus  the  shear  is  also  important  In  the  water  motion  vorti- 
city  constraints. 

The  dynamics  governing  sufficiently  low  (subinertial ) frequency  motions 
in  the  ocean  reduces  to  potential  vorticity  conservation.  A class  of  wave- 
like motions  known  as  topographic  planetary  waves  can  be  trapped  in  areas 

£ 

of  the  ocean  where  topographic  gradients  are  large,  such  as  over  continental 
slopes  and  scarps,  or  mid-ocean  ridges.  Waves  in  this  class  trapped  over 
continental  shelves  and  slopes  are  usually  known  as  continental  shelf  waves 
(CSW's),  as  named  by  Robinson  (1964),  with  coastal  trapped  (boundaries 
important)  or  bottom  trapped  (stratification  important)  variants.  The 
restoring  force  for  CSW's  arises  from  horizontal  gradients  of  the  "back- 
ground" potential  vorticity  field,  f/H,where  f is  suitably  modified  to 
take  account  of  background  current  shears,  and  H is  the  water  depth.  For 
elaboration  of  these  details,  see  Rhines  (1970).  A text  review  of  the 
subject  is  given  by  LeBlond  and  ftysak  (1978). 

Because  of  their  generic  nature,. CSW  horizontal  velocity  (u,v)  and 
sea  level  displacement  (n)  fields  are  trapped  over  steep  bottom  slopes, 
generally  with  oscillatory  cross-slope  structure  in  those  regions  and 
exponentially  decaying  structure  in  regions  of  gentle  slope.  Because  CSW's 
do  not  depend  on  gravity  for  their  existence,  associated  sea  level  displace- 
ments are  small  [0(10  cm)]  and  may  become  undetectable  away  from  the  steeply 
sloping  zone,  because  the  wave  motion  Is  evanescent  there.  Typical  CSW 
velocities  are  0(30  cm  s_1),  with  wavelengths  of  0(100  - 1000  km)  and 


OFFSHORE  OISTANCE  (KM) 


Two-month  summer  mean  longshore  Gulf  Stream  isotachs  off  Cape 
Fear,  after  Richardson  et  at.  (1969).  The  mean  surface  axis 
approximately  overlies  the  550  m isobath. 


6 


"normal"  phase  propagation  cum  sole  (i.e..  In  the  Kelvin  wave  sense)  at 
speeds  of  0(100  km  day"1)*  An  opposing  background  current,  such  as  the 
Gulf  Stream  in  the  SAB  upper  slope  case,  can  strongly  modify  these  charac- 
teristics (Brooks  and  Mooers,  1977a),  and,  under  certain  conditions,  reverse 
the  direction  of  wave  propagation  and  lead  to  wave  Instabilities  (Niiler 
and  Mysak,  1971). 

To  measure  CSW  properties  on  a continental  margin,  one  would  ideally 
employ  current-measuring  techniques  near  the  shelf  edge,  since  the  horizon- 
tal velocity  field  is  the  strongest  CSW  "signal."  As  a relatively  cow  cost 
precursor  to  such  a project,  however,  one  can  take  advantage  of  the  CSW  sea 
level  oscillations  which  exist  in  conjunction  with  the  velocity  field, 
because  they  can  be  detected  in  coastal  tide  gage  records  If  the  recording 
stations  are  sufficiently  close  to  the  "waveguide"  (i.e.,  sufficiently  close 
to  the  steep  bottom  slope,  where  the  wave  is  trapped.  Such  observations 
have  been  made  in  numerous  places  (see  Mysak  (1979)  for  a review)  and  clearly 
indicate  the  importance  of  CSW  motions  over  continental  margins  [eastern 
Australia,  Hamon  (1962);  western  United  States,  Mooers  and  Smith  (1968); 

North  Carolina,  Mysak  and  Hamon  (1969),  and  Brooks  (1978);  Florida,  Brooks 
and  Mooers  (1977a)].  Tide  gages  separated  from  the  shelf  edge  by  wide,  flat 
shelves  are  not  as  useful  as  CSW  detectors  (Brooks  and  Mooers,  1977a),  and  in 
some  cases,  may  be  isolated  from  motions  at  the  shelf  edge  (Lee  and  Brooks, 
1979).  In  shallow  water,  near  shore  wind-driven  motions  may  completely 
swamp  CSW  signatures  In  tide  gage  records.  The  lower  bound  on  detectable 
CSW  sea  level  amplitudes  appears  to  be  several  cm;  thus,  to  be  an  effective 
detector,  the  tide  gage  must  be  within  about  one  e-folding  scale  of  the 
shelf  edge.  For  CSW?s  this  scale  is  X/2ir  where  X is  the  wavelength. 

From  the  regional  tide  gage  studies  in  the  SAB  just  mentioned. 
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indications  of  southward  traveling  CSW's  have  been  found  between  Cape 
Hatteras  and  Southport,  North  Carolina,. and  along  the  southern  half  of  the 
east  Florida  shelf.  However,  the  signal  "fades  out"  of  tide  gage  records 
along  the  relatively  wide  and  shallow  Georgia  shelf  which  joins  the  two 
regions  (Lee  and  Brooks,  1979).  In  the  MAB,  the  evidence  seems  quite  con- 
vincing that  CSW's  are  forced  by  longshore  winds  in  the  New  York  Bight  area, 
and  then  propagate  southward  as  free  waves  at  least  as  far  as  Kiptopeake 
Beach  (Wang,  1979). 

The  principal  objective  of  the  present  study  is  to  use  tide  gage 
records  to  see  if  CSW's  can  propagate  "around  the  corner"  at  Cape  Hatteras, 
thereby  linking  the  Bights.  Such  a link  would  provide  a potentially 
important  mechanism  by  means  of  which  storm  or  otherwise  generated  vorticity 
anomalies  in  the  MAB  could  be  introduced  into  the  SAB.  An  additional 
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objective  is  to  examine  the  relationship  between  the  tide  gage  sea  level 
records  and  the  coastal  meteorological  forcing  fields  of  wind  stress  and 
atmospheric  pressure.  The  remainder  of  the  body  of  the  report  discusses: 

(a)  the  data  set  and  the  methods  of  analysis  used,  b)  summary  results  of 
the  data  analysis,  and  c)  implications  concerning  the  importance  of  CSW's 
in  the  SAB  and  MAB  in  light  of  some  simple  model  calculations.  The  support- 


ing  data  set  from  which  the.  summary  figures  were  extracted  are  contained 


in  appendices. 
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2.  Data  Sources  and  Processing  Techniques 


I 


The  raw  data  set  used  for  analysis  in  this  study  was  extracted* from 
the  NOAA  archives  (hourly  coastal  tide  gage  observations.  National  Ocean 
Survey,  Rockville,  Maryland;  three-hourly  coastal  meteorological  observa- 
tions, National  Climatic  Center,  Asheville,  North  Carolina).  The  observa- 
tion stations  were  chosen  to  provide  reasonable  symmetry  with  respect  to 
Cape  Hatteras,  the  "vertex"  station  of  the  array.  An  additional  goal  was 
to  obtain  data  subsets  of  about  four  months  duration  for  winter  and  summer 
meteorological  conditions,  to  permit  a seasonal  comparison  of  meteorologi- 
cal forcing  characteristics  and  effectiveness.  The  raw  subsets  finally 
chosen  as  the  best  compromise  between  the  goals  and  the  available  data  span 
the  periods  18  Apr-21  Aug  74  ("summer")  and  1 Nov  74-15  Mar  75  ("winter"), 
with  start  times  of  all  data  records  of  0100  local  time.  These  choices 
avoided  all  data  gaps,  except  for  the  winter  Garden  City,  South  Carolina 
tide  gage  record,  where  gaps  of  19,  6,  3 and  14  days  occurred  at  random 
Intervals.  The  gaps  were  filled  with  the  mean  value  of  the  record  over 
suitably-determined  intervals  before  and  after  the  gap.  The  resulting  time 
series  preserves  the  phase  coherence  of  "good"  portions,  but  suffers 
spectral  blurring  in  proportion  to  the  amount  of  data  lost,  which  Is  30. 7% 
for  Garden  City. 

The  raw  meteorological  data  used  in  this  study  consisted  of  wind  speed 
and  direction  and  atmospheric  pressure  at  observation  height.  The  wind 
velocity  vectors  were  converted  to  wind  stress  vectors  (t)  , using  the 
formula 

T - -pa  CB|ll|u 

where  p.  Is  the  air  density  at  sea  level,  the  drag  coefficient 

a 
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Cq  ■ 1.5  x 10"*  is  assumed,  u Is  the  wind  velocity  vector,  and  the  minus 
sign  results  In  the  positive  wind  stress  vector  sense  corresponding  with 
the  direction  toward  which  the  wind  blows.  The  wind  stress  vectors  were 


then  resolved  into  offshore  (t  ) and  longshore  (x  ) components  in  rotated 

a y 


coordinate  systems,  with  the  longshore  direction  parallel  to  the  alignment 
of  the  200  m isobath  immediately  offshore  of  the  meteorological  station. 
This  rotation  scheme,  while  in  general  resulting  in  different  rotation 
angles  for  each  station,  emphasizes  the  expected  importance  of  the  conti- 
nental slope  in  the  relation  between  the  wind  and  sea  level  fields. 

The  raw  data  were  then  low-pass  filtered1  to  remove  diurnal  and  semi- 
diurnal tidal  and  inertial  (the  Cape  Hatteras  inertial  period  is  about  21 
hours)  variance,  and  the  resulting  "40HRLP"  time  series  were  decimated  to  a 
common  six-hour  interval,  with  start  and  stop  times  of  0100  and  2200  local 
time,  respectively.  The  low-passed  sea  level  time  series  were  "barometri- 
cally adjusted"  to  remove  the  static  sea  level  height  changes  associated 
with  barometric  pressure  changes;  this  was  accomplished  by  adding  1.01  cm 
to  the  sea  level  height  at  each  station  for  each  mb  of  pressure  increase 
relative  to  the  mean  pressure  at  the  nearest  meteorological  station.  The 
filtered,  "working"  time  series  span  the  periods  22  Apr-17  Aug  74  (summer) 
and  5 Nov  74-11  Mar  75  (winter),  with  record  lengths  of  118  (127)  days  in 
summer  (winter).  Station  names  and  abbreviations  and  coordinate  rotation 
angles  are  given  in  Table  1,  and  station  locations  are  shown  in  Fig.  1. 


1 A 40-hour  quarter-power  (-6  db)  point  low-pass  Lanczos  filter  kernel 
with  33  (97)  weights  was  used  to  filter  the  3-hourly  meteorological  (hourly 
sea  level)  data,  resulting  in  a data  loss  of  4 days  at  the  beginning  and 


end  of  all  time  series.  Filter  attenuation  at  diurnal  and  higher  frequen- 
cies is  >106  . The  standard  filter  programs  are  called  "METFIL"  and 
"SLVFIL."  resDectivel v.  and  are  available  from  the  author. 
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Table  1.  Meterologlcal  and  sea  level  stations  for  the  sunmer  (S,  22  Apr- 
17  Aug  74,  118  days)  and  winter  (W,  5 Nov  74-11  Mar  75)  40-hour 
low  pass  filtered  data  sets.  The  meteorological  station  coordi- 
nate rotation  angles  are  with  respect  to  north,  positive  clock- 
wise. Station  locations  are  shown  In  Fig.  1. 


Station  name 


A.  Meteorological  Stations 


Abbreviation 


Coordinate 
rotation  (degree) 


1 

i 1 

Norfolk,  VA 

NOR 

10 

! 

\ ■ 

E 1 

Cape  Hatteras,  NC 

HAT 

56 

li 

Wilmington,  NC 

WIL 

56 

1 1 

r> 

1 1 

Charleston,  SC 

CHS 

56 

B.  Sea  Level  Stations 

• 

Station  name 

Abbreviation  and 
season  (S/W)  used 

Nearest  meteorological 
station  for  comparison 

r « 

Atlantic,  City,  NJ 

ACY  (W) 

NOR 

Cape  May,  NJ 

MAY  (S) 

NOR 

Kiptopeake  Beach,  VA 

KIP  (S  and  W) 

NOR  ' 

Avon,  NC 

AVN  (S  and  W) 

HAT 

Beaufort,  NC 

BFT  (S) 

HAT 

Ocean  City  Beach,  NC 

OCY  (W) 

WIL 

Wilmington,  NC 

WIL  (S  and  W) 

WIL 

. 

Myrtle  Beach,  SC 

MYR  (S) 

WIL 

1 ' 

Garden  City,  SC 

GCY  (W) 

WIL 

Charleston,  SC 

CHS  (S  and  W) 

CHS 

* 30.7%  gappy 
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Spectrum  estimates  were  computed  from  the  working  data  set  by  Fourier 
transforming  cosine-tapered  correlation  functions  (Blackman  and  Tukey,  1958). 

K I I 

The  resulting  estimates  covered  the  frequency  range  from  0 to  0.5  cycles  per 
day  (CPD)  with  an  effective  bandwidth  of  0.033  CPO  and  16  (17)  degrees  of 
freedom  for  the  summer  (winter)  case.  Coherence  squared  (y2)  and  phase 

- 

L.  1 

(e)  functions  between  two  time  series  were  computed  from  the  (complex) 

' cross  spectrum  squared  normalized  by  the  product  of  the  autospectra.  For 

* brevity,  y2  will  be  referred  to  as  "coherence."  The  (coherence,  phase) 

vector  Is  a measure  of  linear  dependence  between  two  time  series,  expressed 
as  a complex  function  of  frequency. 

A frequency-dependent  multiple  regression  technique  (Groves  and  Hannan, 
1968)  was  used  to  investigate  meteorological  forcing  of  sea  level.  The 
"Input"  variables  were  the  two  wind  stress  components  and  atmospheric 
pressure.  Multiple  coherence  between  different  combinations  of  Inputs  and 
I • sea  level,  and  complex  partial  transfer  functions  between  Individual  inputs 

and  sea  level  were  computed.  This  technique  has  the  advantage  of  accounting 

I for  possible  input  variable  correlations,  so  that,  for  example,  the  trans- 

fer function  between  atmospheric  pressure  and  sea  level  can  be  corrected 
for  coherence  between  wind  stress  and  pressure. 
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3.  Summary  Analysis  Results 

The  spectrum  and  regression  analysis  results  In  condensed  format  are 
presented  in  this  section  as  contoured  fields  in  the  w-y  plane,  where 
u)  is  frequency  and  y Is  the  distance  from  the  Avon  sea  level  station 
measured  along  the  200  m isobath. 

3.1  Sea  level  interstation  analysis 

The  coherence  and  phase  of  adjusted  sea  level  at  all  stations  relative 
to  AVN  Is  shown  in  Figs.  4-7  for  winter  and  in  Figs.  8-11  for  summer.  The 
95%  null  hypothesis  level  for  coherence  is  about  0.4,  and  the  95%  confidence 
interval  for  phase  is  about  ±(3,  12,  18)  degrees  for  coherences  of  (0.9, 

0.8,  0.7),  respectively.  The  sign  convention  for  the  phase  is  such  that 
positive  (negative)  values  indicate  fluctuations  leading  (lagging)  with 
respect  to  AVN. 

The  adjusted  sea  level  coherence  in  both  seasons  (Figs.  4 and  8)  shows 
selective  interstation  coupling  over  both  Bights  in  period  bands  of  approxi- 
mately 2.5  to  3.5  and  _>  10  days.  The  relative  coherence  maxima  in  these 
bands  can  be  traced  over  most  of  the  extent  of  the  array.  In  the  winter 
2.5  to  3.5  day  band,  for  which  the  coherence  decreases  rapidly  north  of  KIP. 
The  coherences  are  generally  lower  for  the  unadjusted  sea  level  (Figs.  6 and 
10),  except  for  periods  greater  than  about  one  week.  An  additional  band  of 
high  coherence  occurs  in  the  MAB  unadjusted  sea  level  at  a period  of  about 
six  days,  corresponding  to  the  mean  cold  front  interval  for  the  Eastern 
United  States  coastal  margin  (Mooers,  et  at.,  1976).  This  effect  Is  most 
noticeable  in  the  winter  and  In  the  MAB,  as  one  would  expect.  The  baro- 
metric adjustment  procedure  effectively  removed  the  six-day  sea  level 
response,  as  evidenced  by  a comparison  of  Figs.  4 and  6.  In  the  summer. 


n_i . 
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Contours  of  winter  coherence  (x  10)  relative  to  AVN  for  unad 
justed  sea  level  ( cf . Fig.  4). 


Contours  of  winter  phase  (r  10)  relative  to  AVN  for  unadjusted 
sea  level.  The  right  panel  was  contoured  without  the  WIL 
phase  (of.  Fig.  5). 


Contours  of  summer  phase  (t  10)  relative  to  AVN  for  adjusted  sea 
level.  The  right  panel  was  contoured  without  the  WIL  phase 
[of.  Fig.  5). 


<SAva)  aoid3d 


ca  do)  AON3no3dd 
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when  cold  front  passages  over  the  array  area  are  Infrequent,  the  adjustment 
process  essentially  resulted  In  a broad^banded  coherence  Increase  In  both 
Bights  (Figs,  8 and  10). 

The  winter  and  summer  phase  results  for  both  the  adjusted  (Figs.  5 and 
9)  and  unadjusted  (Figs.  7 and  11)  sea  level  are  consistent  with  disturbances 
traveling  southwestward  In  the  high  coherence  bands  over  the  entire  array. 

The  Indicated  propagation  speeds  In  the  2.5  to  3.5  day  band  range  from  a 
low  of  several  hundred  km  d"1  near  AVN  to  essentially  Infinite  (within  the 
contouring  resolution  of  the  figure)  in  the  southern  half  of  the  SAB.  Indi- 
cated speeds  In  the  > 10-day  band  are  of  the  order  of  several  hundred  km  d-1, 
but  the  data  resolution  Is  Insufficient  to  justify  an  Interbight  or  sea- 
sonal comparison. 

The  Wilmington  sea  level  fluctuations  have  anomalously  large  phase 
lags,  as  evidenced  by  the  tendency  toward  vertical  contours  near  WIL  in  the 
center  panels  of  the  phase  figures  (Figs.  5 and  9).  The  anomaly  Is  45-60 
degrees,  or  about  3/8  to  1/2  day  In  the  2.5  to  3.5  day  band,  and  may  be  due 
to  a river  delay  (the  gage  Is  located  about  30  km  upstream  from  the  Cape 
Fear  river  mouth).  However,  lagging  phases  In  the  southwestward  direction 
are  confirmed  without  using  the  WIL  data  In  the  contours  (right-hand  phase 
panels  In  Figs.  5,  7,  9,  and  11). 

The  above  analysis  suggests  that  the  MAB  and  SAB  are  selectively  coupled 
In  the  2.5  to  3.5  and  10-day  period  bands.  However,  the  2.5  to  3.5  day 
band  has  a longer  coherence  scale  In  the  SAB  than  In  the  MAB.  It  Is  also 
apparent  from  the  coherence  figures  (a  clear  example  being  Fig.  4)  that  the 
coherences  relative  to  AVN  are  higher  and  better  organized  (smoother  contours) 
In  the  MAB  than  In  the  SAB.  Part  of  the  relative  organization  In  the  MAB 
may  simply  be  a result  of  fewer  stations  to  contour  there.  However,  AVN 
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may  not  be  well  suited  as  a vertex  station,  because  It  biases  the  analysis 
toward  the  MAB,  since  the  AVN  tide  gage  is  located  about  15  km  north  of 
Cape  Hatteras  and  Is  partially  sheltered  from  the  SAB  by  Diamond  Shoals. 

A better  location  for  a vertex  station  would  probably  be  the  Diamond  Shoals 
light  tower. 

Plotting  the  adjusted  sea  level  coherences  between  adjacent  stations 
shows  that  the  sea  level  fluctuations  differ  substantially  between  the 
Bights,  especially  in  the  winter  (Fig.  12,  winter;  Fig.  13,  summer).  Gen- 
erally the  coherences  are  highest  at  long  periods  in  the  MAB  and  lowest  at 
long  periods  In  the  SAB.  The  SAB  coherences  generally  show  repeatable 
structure,  with  resolved  peaks  occurring  In  both  winter  and  summer.  The  CHS 
record  has  the  highest  coherences  in  both  seasons,  possibly  because  the  sea 
level  amplitudes  at  CHS  are  relatively  large  (see  Appendix  A).  The  MAB 
coherences  show  simpler  structure,  the  principal  feature  north  of  KIP  being 
a sudden  drop  in  coherence  for  periods  shorter  than  about  three  days. 

The  adjacent-station  coherence  analysis  suggests  that  a more  appropri- 
ate SAB  "vertex"  station  is  OCY  in  the  winter  and  BFT  In  the  summer  {of. 

Fig.  1),  thus  avoiding  the  MAB  bias  of  the  AVN  tide  gage.  The  results  of 
this  calculation  are  shown  in  Fig.  14  (winter)  and  Fig  15  (sumner).  In  which 
the  SAB  coherences  are  with  respect  to  the  SAB  vertex  station  and  the  MAB 
coherences  with  AVN  are  reproduced  from  Figs.  4 and  8.  The  SAB  sea  level 
field,  when  examined  Independently  of  the  MAB  field,  appears  to  be  organized 
In  three  bands,  the  periods  of  which  vary  considerably  with  longshelf  posi- 
tion and  with  season.  The  longest  SAB  coherence  scales  (of  order  500  km) 
occur  for  periods  of  about  2.5  and  3.0  days  In  the  winter  and  sumner, 
respectively.  In  contrast,  the  highest  MAB  coherences  occur  for  very  long 
periods  (=  several  weeks)  motions,  especially  In  the  winter. 
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Interstation  winter  adjusted  sea  level  coherence  by  pairs  In  each 
Bight.  The  95%  significance  (null  hypothesis)  level  is  shown  by 
the  horizontal  arrow. 
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Interstation  summer  adjusted  sea  level  coherence  by  pairs  In  each 
Bight.  The  95%  significance  (null  hypothesis)  level  Is  shown  by 
the  horizontal  arrow. 
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Fig.  15  Contours  of  summer  adjusted  sea  level  coherence  (x  10),  relative 
to  AVN  for  MAB  stations  and  relative  to  OCY  for  SAB  stations.. 
Dashed  lines  in  Figs.  15  and  21  to  25  trace  coherence  ridges. 
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The  SAB  sea  level  array  phase  relative  to  the  SAB  vertex  station  (not 
shown)  Indicated  a WIL  river  lag  of  45-60  degrees  in  the  coherent  bands 
with  periods  shorter  than  about  one  week,  as  before.  The  other  SAB  stations 
were  not  significantly  out  of  phase  with  the  SAB  vertex  station  (for  a 
coherence  of  0.6,  the  95%  confidence  interval  on  phase  Is  about  ± 20  degrees), 
although  in  most  cases  small  (of  order  10  degrees)  southwestward  lags  were 
Indicated  In  both  seasons.  Reexamining  Figs.  5 and  9,  one  sees  that  most 
of  the  SAB  phase  lag  relative  to  AVN  occurs  between  AVN  and  the  SAB  vertex 
station,  and  that  the  remaining  SAB  stations  (WIL  excepted  because  of  the 
river  lag)  are  nearly  in-phase  in  the  coherent  bands. 

3.1.1  Discussion.  The  low  frequency  sea  level  fluctuations  in  both 
Bights  are  organized  in  several  distinct  period  bands,  with  longshelf 
coherence  scales  of  at  least  500  km  in  the  SAB  and  400  km  in  the  MAB.  How- 
ever, there  are  seasonally  dependent  interbight  differences  between  the 
periods  of  the  coherent  bands.  The  clearest  picture  emerges  when  each 
bight  Is  considered  separately,  suggesting  Important  differences  in  the  pro- 
cesses affecting  the  SAB  and  MAB  coastal  sea  level  and  therefore  (by  Infer- 
ence) affecting  the  water  motions  over  the  continental  margin.  Three 
candidate  explanations  for  the  differences  between  the  SAB  and  MAB  are: 
a)  differences  In  meteorological  forcing  characteristics,  b)  differences  in 
the  continental  margin  topography,  and  c)  the  presence  of  the  Gulf  Stream 
in  the  SAB.  The  first  possibility  is  considered  in  Section  3.2  and  the 
latter  two  are  addressed  in  Section  4.0. 

The  coherence  structure  noted  here  is  consistent  with  that  found  in 
several  previous  sea  level  investigations  In  the  MAB-SAB  area.  In  a study 
of  BFT  and  WIL  data.  Brooks  (1978)  found  three  distinct  peaks  in  the  low 
frequency  coherence  between  the  barometrically  adjusted  sea  level  records. 


which  he  interpreted  as  near-zero  group  speed  CSW  responses  to  wind  forcing. 
In  an  MAB  sea  level  study  extending  from  AVN  to  Nantucket  Shoals  (about  400 
km  northeast  of  MAY  in  Fig.  1,  Wang  (1979)  found  only  very  long  period  co- 


herent motions  north  of  KIP,  with  evidence  of  short  period  motions  appear- 


ing between  AVN  and  KIP.  Thus  it  seems  that  subtidal  coastal  sea  level 


motions  with  periods  shorter  than  about  a week  are  a feature  of  the  SAB  but 


are  evanescent  north  of  AVN. 


The  most  prominent  feature  of  the  longshelf  sea  level  phase  analysis 


is  the  WIL  tide  gage  river  lag  of  45-60  degrees,  equivalent  to  a time  lag 


of  about  3/8  day  in  the  three-day  period  band.  The  other  station  phases 


generally  show  small  lags  in  a southwestward  direction,  but  they  are  only 


questionably  different  from  zero  at  the  95%  confidence  level.  Thus  the 


evidence  for  longshelf  propagation  of  disturbances  is  murky  at  best  and 


is  not  easily  explained  in  terms  of  a single  propagating  wave.  It  is  now 


also  apparent  that  the  BFT  to  WIL  phase  lag  reported  by  Brooks  (1978)  may  be 


partly  accounted  for  as  Cape  Fear  river  lag,  a possibility  raised  in  that 


paper,  and  that  in  fact  the  sea  level  fluctuations  at  BFT  and  at  the  river 


mouth  are  nearly  in-phase.  In  an  earlier  study  of  BFT  and  Southport,  NC 


(at  the  Cape  Fear  river  mouth)  tide  gage  data,  Mysak  and  Hamon  (1969)  found 


maximum  low  frequency  lags  of  about  10  degrees,  consistent  with  the  south- 


westward  traveling  disturbances,  but  implying  propagation  speeds  consider- 


ably larger  than  expected  from  a free  CSW  model.  Their  10  degree  lag  w« 


barely  significant  at  th£  95%  level,  because  a full  year's  data  set  was  used 


for  the  analysis.  In  view  of  the  consistency  of  the  southwestward  but  small 


lags  over  the  present  array,  the  statistical  significance  test  based  on  only 
one  "realization  of  the  process"  (i.e.,  on  only  one  tide  gage)  may  be 


unduly  conservative.  If  one  adopts  this  view,  it  can  be  argued  that  the 
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phase  results  Indicate  southwestward  traveling  disturbances  over  most  of 
the  array  in  the  coherent  bands.  However,  the  very  small  lags  (tens  of 
degrees)  involved  lead  to  implied  wavelengths  many  times  longer  than  the 
longshelf  coherence  scales,  whence  the  concept  of  a traveling  monochromatic 
wavelike  disturbance  becomes  untenable.  A more  attractive  alternate  inter- 
pretation  may  be  that  locally  forced  disturbances  overwhelm  remotely  forced 
disturbances,  as  suggested  by  Mysak  and  Hamon  (1969).  A related  possibility 

| 

is  that  wave  group  properties  may  be  more  appropriate  than  phase  properties 
when  trying  to  explain  the  coastal  sea  level  fluctuation  field. 

3.2  Meteorological  forcing 

In  this  section,  a frequency-domain  linear  regression  technique  is  used 
to  investigate  the  combined  and  individual  effects  of  atmospheric  pressure 
and  wind  stress  on  sea  level  height.  Since  the  pressure  and  wind  stress 
forcing  contributions  can  be  expected  to  be  correlated  to  some  degree,  it  is 
useful  to  consider  their  collective  influence  using  a multiple  input  stati- 
cal model  which  corrects  individual  transfer  functions  for  the  coherent  con- 
tributions of  other  inputs.  Thus,  for  example,  the  transfer  function  between 
pressure  and  sea  level  fluctuations  can  be  corrected  in  each  frequency  band 
to  remove  the  effects  of  coherent  wind  stress  (Bendat  and  Piersol,  1968; 

Groves  and  Hannan,  1968). 

Meteorological  (MET)  time  series  of  pressure  (p)  , cross  shelf  wind 
stress  (xu)  and  longshelf  wind  stress  (xy)  at  NOR,  HAT,  and  WIL  were 

used  as  forcing  inputs  for  five  winter  and  four  summer  sea  level  stations  . 

' 

spanning  both  Bights  (see  Table  l.B  for  correspondence  between  meteorologi- 
cal and  sea  level  stations).  The  results  are  contoured  in  w-y  space,  and 
supporting  calculations  are  shown  in  Appendix  C.  For  each  season,  the 
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following  quantities  are  discussed:  a)  the  multiple  coherence 

y2_  ci  ) which  is  a measure  of  the  simultaneous  linear  dependence 

P*  u*  V* 

of  the  sea  level  (SL)  fluctuations  on  all  three  input  variables;  b)  the 
pressure  and  wind  stress  component  partial  transfer  function  moduli 

^HP*SL;tu,tv^  * ^Htu-SL;p,tv^  * ^tv*SL;p,tu)  , where  the  transfer  func- 
tion between  the  two  dotted  subscript  variables  has  been  corrected  for 
remaining  ones;  and  c)  the  coherent  fraction  of  longshelf  wind  stress  com- 
pared to  total  wind  stress  [TV/(TU+TV) 1 » which  is  a measure  of  the 

importance  of  wind  forced  coastal  Ekman  flux  divergence.  The  interstation 
coherence  ridges  identified  in  Figs.  14  and  15  have  been  transferred  to 
the  figures  in  this  section  for  comparison. 

The  WIL  transfer  functions  all  had  anomalously  large  moduli,  indicative 
of  shallow  water  or  closed  basin  responses  atypical  in  magnitude  from 
coastal  responses.  For  this  reason,  the  WIL  moduli  were  not  used  in  the 
contours,  although  frequencies  of  relative  maximum  and  minimum  values  often 
were  identifiable  with  those  at  coastal  stations.  The  WIL  sea  level  record, 
while  evidently  biased  by  river  lag  and  shallow  basin  effects,  still  contains 
a useful  "signature"  of  coastal  processes. 

3.2.1  The  winter  case.  The  multiple  coherence  (Fig.  16)  indicates 

•k 

significant  coupling  between  MET  and  SL  fluctuations  over  a broad  area  of 
w-y  space  in  both  Bights.  The  area  includes  all  the  coherent  interstation 
period  bands  identified  earlier  (the  dashed  lines)  except  the  short  period 
MAB  band,  for  which  the  coherence  falls  well  below  the  significance  level. 

The  95%  significance  level  for  y*  _ .Cl 

P»Tu,Tv  3L 


is  0.45. 


This  supports  the  notion  raised  earlier  that  short  period  motions  can  be 
forced  In  the  SAB  but  are  evanescent  in  the  MAB. 

The  non-static  residual  partial  pressure  transfer  function  is  shown  In 
Fig.  17.  The  residual  was  obtained  by  subtracting  1.01  cm/mb  (the  static 
"barometer  factor")  from  each  transfer  function  estimate.  Thus  if  the  sea 
level  response  to  atmospheric  pressure  were  entirely  static,  the  residual 
would  be  zero;  conversely,  large  departures  from  zero  of  either  sign  in 
Fig.  17  indicate  areas  of  large  dynamical  sea  level  responses.  Comparing 
Figs.  16  and  17,  one  sees  considerable  alignment  of  the  coherent  sea  level 
bands  and  the  coherent  MET’SL  areas  with  areas  of  dynamic  pressure  trans- 
fer function,  indicating  that  the  atmospheric  pressure  makes  an  important 
contribution  to  the  winter  forcing.  For  periods  shorter  than  a week,  the 
pressure  forcing  is  confined  almost  entirely  to  the  SAB.  It  is  also  note- 
worthy that  the  highest  MET*SL  coherences  occur  in  areas  of  positive  (super- 
barometric)  pressure  transfer  function  residual. 

The  cross  shelf  and  longshelf  wind  stress  transfer  functions  (Figs.  18 
and  19,  respectively)  show  that  effective  wind  stress  forcing  of  sea  level 
occurs  in  both  Bights  at  periods  longer  than  about  one  week,  contrasting 
with  the  pressure  case  in  Fig.  17.  The  wind  stress  components  are  roughly 
equally  important  except  between  OCY  and  AVN,  where  the  longshelf  component 
dominates  for  periods  longer  than  a week.  The  short  and  intermediate  period 
coherent  sea  level  bands  (dashed  lines)  do  not  seem  to  be  identifiable  with 
features  of  the  wind  stress  forcing  field.  For  long  periods,  however,  both 
wind  stress  components  contribute  to  the  MET*SL  coherence  in  Fig.  16;  the 
large  ry  transfer  function  values  are  clearly  associated  with  the  sea 
level  coherence  maximum  near  AVN.  It  also  can  be  seen  by  comparing  Fig.  16 
and  the  coherent  longshelf  wind  stress  fraction  in  Fig.  20  that  the  areas 
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of  highest  MET*SL  coherence  are  not  aligned  with  areas  of  highest  values 
of  the  fraction;  l.e.,  the  cross  shelf  wind  stress  plays  an  important  role 
In  the  forcing  process.  A similar  effect  was  noted  in  the  Florida  Straits 
by  Brooks  and  Mooers  (1977a).  Thus  the  coastal  Ekman  divergence  shelf  wave 
forcing  mechanism  described  by  Gill  and  Schumann  (1974)  may  not  be  entirely 
adequate  here,  although  It  Is  clear  from  Fig.  20  that,  overall,  xy  is  more 
I effectively  coupled  to  sea  level  fluctuations  than  is  tu  (the  area-weighted 

fraction  in  Fig.  20  is  about  0.8). 

3.2.2  The  summer  case.  The  suimer  MET*SL  coherence  (Fig.  21)  has 

more  and  larger  areas  of  high  coherence  (>  0.8)  than  the  winter  case.  The 

short  and  intermediate  period  sea  level  bands  (dashed  lines),  while  not  as 
clearly  associated  with  high  MET*SL  coherence  as  they  were  In  the  winter, 
fall  within  the  5.5  contour  (shown  as  a special  case  In  Fig.  21),  and  the 

long  period  band  In  both  Bights  falls  Inside  or  near  the  0.8  contour. 

Comparing  Fig.  21  with  the  residual  pressure  transfer  function  In 
Fig.  22,  one  sees  correspondence  between  high  MET  SL  coherence  and  dynamic 
pressure  forcing  In  the  SAB  for  periods  shorter  than  a week,  but  little 
correspondence  elsewhere  In  the  w-y  space.  Overall,  the  simmer  sea  level 
pressure  response  was  more  nearly  static  than  In  the  winter,  but  with  the 
largest  residuals  occurring  for  short  periods  In  the  SAB  as  in  the  winter. 

The  MAB  coherent  sea  level  fluctuations  appear  to  be  unrelated  to  dynamic 
pressure  fluctuations.  Both  wind  stress  components  were  relatively  more 
effective  In  the  SAB  than  In  than  In  the  MAB  for  periods  longer  than  about 
one  week  (Figs.  23  and  24),  but  the  maximum  magnitudes  of  both  components 
were  about  30%  less  than  In  the  winter.  As  In  the  winter,  the  high  MET*SL 
coherence  areas  are  anti -correlated  with  areas  of  high  coherent  ry  frac- 
tion (Fig.  25),  again  indicating  the  Importance  of  the  cross  shelf  wind 
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Contours  of  summer  multiple  coherence  (*  10)  between  sea  level 
and  the  combined  effects  of  atmospheric  pressure,  cross  shelf 
wind  stress,  and  long  shelf  wind  stress. 
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stress  component. 

3.2.3  Discussion.  There  are  several  seasonally  persistent  charac- 
teristics of  the  meteorologically  forced  sea  level  field  common  to  both 
Bights.  Strong  coupling  of  the  wind  stress  and  sea  level  fluctuations  occurs 
only  for  periods  longer  than  about  one  week,  and  both  the  longshelf  and  the 
cross  shelf  wind  stress  components  are  important  forcing  contributors  in 
the  coherent  sea  level  bands.  The  dynamical  pressure  forcing  contribution 
is  important  for  short  period  motions,  especially  noticeable  in  the  SAB,,  and 
more  prominently  so  in  the  winter  than  in  the  summer.  On  the  other  hand,  the 
effective  wind  stress  forcing  in  the  summer  is  largely  confined  to  the  SAB, 
whereas  in  the  winter  it  occurs  in  both  Bights.  This  result  probably 
reflects  seasonal  differences  in  weather  patterns,  with  cold  front  winds 
influencing  both  Bights  in  the  winter  and  the  more  gentle  and  steady  south- 
east trade  winds  primarily  affecting  the  SAB  in  the  summer. 

For  several -day  periods,  the  coastal  sea  level  coherence  pattern  seems 
to  be  associated  with  dynamical  sea  level  responses  to  atmospheric  pressure 
fluctuations,  but  only  in  the  SAB.  Strong  dynamical  sea  level  responses  to 
atmospheric  pressure  were  also  found  off  Australia  by  Hamon  (1966),  off 
North  Carolina  by  Mysak  and  Hamon  (1969)  and  off  Florida  by  Brooks  and  Mooers 
(1977a),  suggesting  the  general  importance  of  pressure  forcing.  Adams 
and  Buchwald  (1968)  have  argued  that  the  correlation  between  sea  level  and 
pressure  observed  by  Mysak  and  Hamon  is  an  arti factual  consequence  of  corre- 
lation between  the  pressure  and  the  longshelf  wind  stress,  with  the  latter 
providing  the  "real"  vorticity  forcing  via  coastal  Ekman  flux  divergences. 

The  pressure  transfer  functions  presented  herein,  however,  have  been  cor- 
rected for  the  linear  correlation  with  wind  stress,  and  the  non-static 
responses  are  still  clearly  evident.  It  is  furthermore  evident  that  the 


pressure  and  wind  stress  forcing  contributions  generally  occupy  different 
portions  of  the  frequency  domain;  i.e.,  the  pressure  effect  is  relatively 
more  important  at  short  periods,  and  the  wind  stress  effect  is  more  impor- 
tant at  long  periods. 

The  argument  that  pressure  fluctuations  are  ineffective  vorticity 
generators  over  a shelf  is  based  on  the  fact  that  the  free  surface  diver- 
gence term  in  the  barotropic  long  wave  equations  is  scaled  by  (L/Rq)2  , 
where  L is  the  cross  shelf  length  scale  and  Rq  is  the  Rossby  deforma- 
tion radius.  For  typical  shelves,  this  scale  factor  is  0(5  x 10"2);  thus, 
in  the  absence  of  resonances,  the  above  argument  follows.  Adams  and  Buch- 
wald  point  out  that  a resonant  response  with  the  pressure  disturbance 
traveling  in  phase-locked  coincidence  with  the  sea  level  disturbance  is  not 
likely  to  be  important,  because  of  the  long  (1000's  of  km)  coastal  distances 
over  which  the  process  must  act  to  be  effective.  How  then  are  the  dynamic 
sea  level  responses  to  pressure  fluctuations  explained? 

The  phase-locked  resonance  mechanism  was  suggested  by  Robinson  (1964). 
An  alternate  type  of  forcing  that  may  be  important  for  both  pressure  and 
wind  stress,  suggested  by  Buchwald  and  Adams  (1968),  may  occur  at  frequen- 
cies for  which  the  group  velocity  of  CSW's  is  small  or  zero.  Since  a CSW 
packet  is  stationary  or  nearly  so  in  the  vicinity  of  2ero  group  velocity 
frequencies  (turning  points  of  the  CSW  dispersion  curves),  the  energy 
supplied  locally  to  such  a packet  of  waves  must  remain  in  its  neighborhood 
of  introduction,  allowing  the  local  packet  amplitude  to  grow  as  long  as  the 
forcing  energy  source  is  available.  Phase  locking  of  the  forcing  and  shelf 
water  disturbances  over  large  coastal  distances  is  not  required  for  this 
"group  resonance"  to  occur.  The  group  resonance  mechanism  also  provides  a 
potential  explanation  for  the  selectivity  (i.e.,  bandedr.ess)  of  the  sea 
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level  coherence  structure.  Several  studies  have  indicated  that  peaks  in 
the  longshelf  coherence  between  sea  level  stations  can  be  associated  with 
the  zero  group  velocity  frequencies  of  CSW  modes  [of.  Cutchin  and  Smith, 
1973;  Brooks,  1978). 

The  wind  stress  and  its  curl  are  forcing  terms  that  may  be  of  compar- 
able magnitude  for  typical  event-like  cold  front  wind  fields  (Brooks,  1978). 
Historically,  the  curl  term  has  been  neglected  compared  to  the  stress  term, 
based  on  the  assumption  that  the  horizontal  scale  of  wind  variations  is 
much  larger  than  the  shelf  width  (see  LeBlond  and  Mysak  for  details  of 
scaling  the  barotropic  vorticity  equation).  However,  it  has  been  pointed 
out  by  Brooks  (1978)  and  others  that  the  horizontal  scale  of  wind  stress 
variation  near  a cold  front  can  be  much  smaller  than  the  distance  between 
adjacent  fronts;  a typical  curl  scale  is  100  km,  the  same  as  the  shelf  width 
scale.  Furthermore,  the  resultant  pulse-like  time  history  of  a sequence  of 
cold  front  wind  stress  curl  events  at  a fixed  location  produces  a forcing 
spectrum  harmonically  rich  in  frequency  and  wavenumber  space.  Thus,  while 
the  fundamental  atmospheric  forcing  time  and  space  scales  may  be  a week  and 
a megameter,  significant  forcing  energy  can  also  be  expected  to  occur  at 
scales  shorter  than  these, especially  for  the  wind  stress  curl.  A field 
test  of  the  importance  of  this  idea  awaits  higher  resolution  wind  stress 
measurements  permitting  a reasonable  calculation  of  the  local  (shelf  scale) 
curl . 

The  longshelf  wind  stress  can  Induce  cross  shelf  water  motions  and 
coastal  sea  level  set-down  or  set-up  in  response  to  surface  Ekman  flux 
divergences  (coastal  or  otherwise  generated),  which  "releases"  vorticity  by 
moving  fluid  columns  over  variable  depth.  However,  for  a straight  coastline 
and  parallel  Isobaths,  a uniform  cross  shelf  wind  stress  can  produce  only 
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a uniform  longshelf  current  geostrophically  balanced  with  a coastal  sea 
level  set  down  or  set  up.  This  situation  would  not  produce  relative  vorti- 
city  anomalies,  assuming  the  longshelf  scales  are  small  enough  that  varia- 
tions in  f (the  "beta-effect")  can  be  neglected,  and  therefore  waves  would 
not  be  excited.  If,  however,  longshelf  coastal  and/or  isobath  variation  is 
significant,  or  if  the  cross  shelf  wind  stress  is  not  uniform,  then  long- 
shelf mass  convergences  and  divergences  can  occur,  resulting  in  cross  shelf 
(vorti city-generating)  water  motions  to  satisfy  mass  conservation  locally. 
Thus  it  seems  that,  while  the  longshelf  wind  stress  mechanism  is  more  direct, 
forcing  by  the  cross  shelf  component  is  likely  to  be  important  in  regions 
of  convergent  or  divergent  isobaths,  or  in  regions  where  isobaths  and  the 
coastline  are  non-parallel,  such  as  in  the  SAB. 

The  results  of  this  section  indicate  that  atmospheric  pressure  and  both 
wind  stress  components  are  important  contributors  to  the  forcing  of  low  fre- 
quency coastal  sea  level  motions.  Their  relative  effectiveness  varies 
according  to  frequency,  season,  and  longshelf  location,  but  it  is  apparent 
that  none  should  be  neglected  a priori  as  insignificant.  It  is  also  clear 
that  the  three-input  forcing  model  presented  does  not  consider  all  possibly 
relevant  contributors,  the  most  notable  absence  being  the  wind  stress  curl. 
The  surprising  degree  of  selectivity  of  the  coupling  between  the  atmosphere 
and  sea  level  fluctuations  suggests  the  importance  of  some  kind  of  resonant 
dynamical  mechanism. 
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4.  Dispersion  of  Barotropic  Shelf  Waves  in  the  Bights 

It  was  pointed  out  in  the  introduction  that  the  continental  slope  and 
shelf  width  change  by  a factor  of  2 or  3 between  Capes  Remain  and  Hatteras. 

A comprehensive  theory  of  coastal  trapped  waves  in  the  SAB  would  therefore 
necessarily  be  three  dimensional  and  would  probably  require  a complex 
numerical  scheme  to  obtain  solutions.  The  presence  of  the  mean  baroclinic 
Gulf  Stream,  which  crosses  bottom  contours  in  the  northern  end  of  the  SAB, 
further  complicates  matters. 

The  general  problem  of  CSW's  propagating  in  a stratified  ocean  with  a 
sheared  mean  current  and  longshelf  variation  of  geometry  has  been  approached 
from  several  simplified  starting  points.  (For  a thorough  review  of  histori- 
cal and  recent  progress,  see  Mysak,  1979).  With  longshelf  variation  neg- 
lected, and  in  the  absence  of  stratification,  continental  shelf  waves 
result  (Robinson,  1964).  When  stratification  is  included,  the  more  general 
coastal  trapped  waves  result  (Wang  and  Mooers,  1976).  If  a vertical  coastal 
wall  is  present,  Kelvin  waves  appear  in  both  cases  which  can  interact  with 
and  modify  the  dispersion  properties  of  the  continental  shelf  or  coastal 
trapped  waves  (Wang  and  Mooers,  1976).  For  the  case  of  small  amplitude 
longshelf  topographic  and  coastline  variations,  barotropic  shelf  waves 
propagate  with  the  local  phase  speed  and  adjust  their  cross  shelf  modal 
Structure  to  agree  with  the  local  eigenfunction  structure  (Allen,  1976). 
During  the  adjustment  process,  scattering  of  energy  into  other  shelf  modes 
can  occur.  If  a sheared  mean  current  flows  parallel  to  isobaths,  the  shear 
and  advective  effects  combined  can  lead  to  stable  shelf  waves  propagating 
in  either  longshelf  direction,  or  to  downstream  unstable  waves  (Niiler  and 
Mysak,  1971;  Brooks  and  Mooers,  1977b).  Much  work  remains  to  be  done  for 
cases  which  Include  large  longshelf  topography  variations  and  baroclinic 
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mean  currents,  particularly  with  regard  to  critical  level  wave-mean  current 
interactions  (of.  McKee,  1979). 

As  a precursor  to  a three-dimensional  formulation,  the  SAB  problem  is 
addressed  in  this  section  with  barotropic  model  in  which  longshelf  varia- 
tions are  assumed  to  occur  over  distances  much  greater  than  the  shelf  width. 
Thus  we  can  expect  slow  variations  in  wave  propagation  speeds  and  eigen- 
functions, as  described  by  Allen  (1976).  The  barotropic  model  discussed  by 
Buchwald  and  Adams  (1968),  modified  to  include  a mean  current,  will  be 
applied  to  bottom  and  mean  current  profiles  representative  of  each  SAB  Cape 
and  the  MAB,  but  the  y (longshelf)  variation  is  not  explicitly  included. 
This  method  is  not  entirely  satisfactory  because  scattering  effects  are 
ignored  and  the  scale  over  which  the  topography  changes  is  too  short  to 
justify  a WKB  approach.  Nonetheless,  the  results  suggest  how  to  proceed. 

4.1  Formulation 

In  the  presence  of  a horizontally-sheared  basic  state  current,  linear 
barotropic  oscillations  of  a fluid  overlying  a sloping  bottom  are  governed 
by  the  potential  vorticity  equation 

<V7>  [(£)  O . 

The  mass  transport  stream  function  has  the  wavelike  form  ^ =tl>(x)e^Ay~a1^  , 
and  the  right-handed  coordinate  system  is  oriented  with  the  x-axis  perpen- 
dicular to  isobaths  (positive  seaward).  Longshelf  (y)  variations  in  the 
mean  current  V(x)  and  the  bottom  topography  h(x)  are  not  explicitly 
included.  The  basic  state  potential  vorticity  is  P(x)  = (f  + Vx ) h “ 1 , 
where  f is  the  constant  Coriolis  parameter.  The  absolute  (i.e.,  relative 
to  the  fixed  coordinate  frame)  wave  phase  speed  is  defined  as  cQ  = <*>/£. 
where  w is  the  wave  frequency  and  i,  is  the  wave  number.  For  a dis- 
cussion of  the  assumptions  and  details  leading  to  this  equation  and  its 


Inherent  limitations,  see  LeBlond  and  Mysak  (1979). 

With  two-point  boundary  conditions  for  , the  governing  equation 
presents  an  eigenvalue  problem  for  the  phase  speed  cQ  , (upstream  propaga- 
tion for  a)  > 0 assumed).  For  downstream  waves  with  0 < cQ  £ 7 , the 

equation  has  singular  (critical)  points  where  cQ  = V , the  spectrum  Is  con- 
tinuous (McKee,  1977),  and  under  certain  circumstances  the  waves  may  be 
unstable  (Nlller  and  Mysak,  1971). 

For  a given  h(x)  , V(x)  , there  will  be  n solutions  (cQ »^)n  For 
each  wavenumber  £n  , with  the  nth  solution  determining  a point  on  the 
dispersion  curve  for  the  nin  mode.  For  example,  the  dispersion  diagram  Is 
given  by  Brooks  (1978)  for  the  Cape  Fear  bottom  profile  (Fig.  2)  and  a 
current  profile  determined  by  vertically  Integrating  the  V(x,z)  section 
In  Fig.  3.  For  general  h , 7 profiles,  the  solutions  are  most  conven- 
iently obtained  numerically. 

The  model  calculations  were  performed  for  each  of  the  bottom  profiles 
shown  In  Fig.  2,  assuming  that  the  mean  current  profile  was  invariant  In 
the  SAB,  that  Its  axis  coincided  with  the  550  m isobath,  and  that  its  trans- 
port was  fixed  at  53  x lO'n^s"1  (reported  by  Richardson  et  al.t  1969). 

The  MAB  case  was  run  with  7 - 0 (no  Gulf  Stream).  Dispersion  curves  for 
the  gravest  mode  In  each  case  are  shown  In  Fig.  26.  The  dimensional  period 
scale  shown  applies  only  at  the  Cape  Fear  latitude  but  the  scale  factor  (the 
Inertial  period)  changes  by  only  about  4 percent  over  the  latitude  range 
considered  (about  3 degrees).  The  zero  group  speed,  lowest  mode  horizontal 
velocity  structure  Is  shown  In  Fig.  27  for  the  three  southernmost  Capes. 

4.2  Discussion  and  relation  to  observations 


For  monotonic  bottom  profiles,  the  dispersion  curve  for  CSW  mode  n 
characteristically  exhibits  at  least  one  zero  group  speed  (zero  slope) 
point  (“no**,, 0)  (Huthnance,  1975).  The  dispersion  curves  In  Fig.  26 


1.5  0.6 t ^DELAWARE  BAY 


Fig.  26  Gravest  shelf  wave  dispersion  curves  for  bottom  and  meai 
current  profiles  typical  of  the  four  Carolina  Capes  and 
the  MAB  (Delaware  Bay).  Zero  group  speed  points  are  cl 
The  dimensional  period  scale  is  strictly  correct  only  a 
Fear. 
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Fig.  27  Gravest  mode  zero  group  speed  velocity  component  eigenfunctions 
for  the  Carolina  Capes.  Locations  of  the  200  m and" 1000  m iso- 
baths and  the  model  mean  Gulf  Stream  axis  arc  shown  by  vertical 
arrows.  Note  that  the  first  mode  structure  over  the  narrow 
Cape  Lookout  shelf  divides  Into  two  coupled  first  mode  struc- 
tures over  the  dual  shelf  off  Cape  Romain;  in  contrast  the 
third  Romain  mode  is  trapped  entirely  over  the  inner  shelf. 
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Illustrate  that  the  location  of  these  points  depends  strongly  on  the  bottom 
profile  and  mean  current  details.  There,  are  at  least  three  competing 
Influences  which  affect  the  overall  shape  of  the  dispersion  curves.  First, 
bottom  slope  Increases  tend  to  Increase  uj  and  i (see  Buchwald  and 
Adams  (1968)  for  a simplified  example).  Second,  the  advection  of  the  mean 
current  In  opposition  to  the  "normal"  direction  of  CSW  propagation  tends  to 
decrease  u>n  for  all  £n  (Brooks  and  Mooers,  1977b).  Third,  the  mean 
shear  of  the  current  aids  or  opposes  the  upstream  tendency  in  regions 
where  the  shear  Is  positive  or  negative,  respectively  (Brooks  and  Mooers, 
1977b).  These  effects  can  be  separated  relatively  easily  in  the  case  of  an 
exponential  bottom  profile,  as  used  by  Buchwald  and  Adams,  but  in  the  pre- 
sence of  more  general  bottom  and  mean  current  profiles,  as  used  here,  they 
may  interact  strongly.  For  example,  the  Cape  Romain  topography  has  inner 
and  outer  zones  of  relatively  steep  slope  separated  by  a zone  of  gentle 
slope.  This  results  in  a splitting  of  the  "waveguide"  for  CSW's  with  suffi- 
ciently short  wavelengths.  The  effect  can  be  seen  in  Fig.  27,  where  the 
Cape  Romain  velocity  structure  resembles  two  coupled  first  mode  waves,  one 
concentrated  on  the  inner  slope  and  one  over  the  outer  slope.  However, 
the  mean  current  in  this  area  is  concentrated  over  the  gentle  slope  region 
between  the  inner  and  outer  waveguides.  Therefore,  the  advective  (kinematic) 
modification  is  weakened,  and  the  Cape  Romain  dispersion  curve  In  Fig.  26 
rises  to  a higher  wn0  than  does  the  corresponding  curve  at  Cape  Fear, 
where  the  waveguide  location  more  closely  coincides  with  the  strong  current 
zone.  At  Cape  Hatteras,  the  very  steep  bottom  slope  and  the  weakened  effect 
of  the  barotroplc  mean  current  (due  to  the  greater  depths  but  constant  trans- 
port) result  in  unQ  , lnQ  being  offscale  in  Fig.  26.  Likewise,  in  the 
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MAB  (represented  by  the  Delaware  Bay  curve  in  Fig.  26),  the  absence 
of  the  current  allows  wnQ  and  Ano  to  increase  to  relatively  large 
values. 

The  foregoing  results  suggest  that  longshelf  energy  dispersion  by  CSW's 
over  the  eastern  United  States  continental  margin  will  be  complex.  For 
example,  consider  the  case  of  a forcing  event  occurring  simultaneously  in 
the  SAB  and  the  MAB.  The  event  can  be  characterized  as  an  energy  density 
distribution  In  u , 1 (dispersion)  space,  which  will  project  into  the 
various  free  modes  of  motion  available  (CSW  s in  this  case)  with  efficien- 
cies dependent  upon  the  eigenfunction  structure  and  other  details  involving 
the  othogonality  of  the  modes  and  critical  level  effects  (of.  Wang  and 
Mooers  (1975)  for  a discussion  of  orthogonality  relations  of  coastal 
trapped  waves  and  McKee  (1979)  for  forced  waves  in  the  presence  of  critical 
levels).  It  is  reasonable  to  expect  the  gravest  mode  to  be  most  strongly 
excited  because  it  has  the  simplest  structure.  The  energy  of  the  mode 
propagates  with  the  group  speed  for  that  mode,  which  is  equal  to  the  slope 
of  its  dispersion  curve.  For  &n  < A , the  energy  flux  Is  upstream  (south- 
westward)  and  for  An  > Ano  , the  energy  flux  is  downstream  (Fig.  26).  For 
the  frequency  and  wavenumber  ranges  in  Fig.  26,  which  correspond  approximately 
with  the  ranges  over  which  the  data  studies  in  Section  3 apply,  a downstream 
energy  flux  can  occur  only  for  periods  of  about  2-4  days,  and  only  in  the 
SAB.  Thus  gravest-mode  backscatterlng  of  CSW  energy  by  topography  or  mean 
current  variations  can  be  significant  only  In  this  period  range.  This  may 
provide  an  explanation  for  the  Intensity  of  short  period  sea  level  motions 
in  the  SAB  compared  to  the  MAB,  and  also  for  the  evanescence  of  short  period 
motions  backscattered  Into  the  MAB  as  long  period  CSW’s  encounter  the 
abrupt  waveguide  change  at  Cape  Hatteras  associated  with  the  Gulf  Stream. 
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The  response  energy  at  a fixed  observation  site  (e.g.,  a tide  gage)  to  the 
simultaneous  event  forcing  would  be  composed  of  the  sum  of  the  energy  Intro- 
duced locally  plus  that  propagated  to  the  observation  site  from  upstream  or 
downstream  forcing  sites.  The  resultant  phase  distribution  along  the 
coast  would  be  determined  in  a complex  way  by  the  relative  importance  at 
each  observation  site  of  the  locally  forced,  forward  propagating,  and  back- 
scattered  energy.  The  fact  that  SAB  tide  gage  investigations  (of.  Section 
3)  have  found  longshelf  phase  lags  in  the  several -day  period  band  too  small 
to  conveniently  explain  by  a monochromatic  disturbance  may  be  evidence  that 
the  locally  forced  contribution  to  the  total  energy  at  a particular  point 
is  dominant.  If  so,  the  dominant  sea  level  and  attendant  circulation  vari- 
ability can  be  expected  to  occur  at  or  near  zero  group  speed  periods  and 
wavelengths.  In  the  SAB,  these  periods  fall  In  the  2-3  day  period  range 
for  the  lowest  CSW  mode,  corresponding  rather  well  to  the  short  period 
coherent  sea  level  bands  In  Section  3.  In  the  MAB,  however,  the  zero  group 
speed  wavelength  Is  very  short  for  the  lowest  mode  — probably  less  than 
the  shelf  width  — and  Is  not  likely  to  be  strongly  energized  by  atmospheric 
disturbances,  even  when  the  relatively  small  wind  stress  curl  scales  (a  few 
hundred  km)  are  considered.  Thus  the  dominant  disturbances  In  the  MAB 
are  likely  to  be  those  with  scales  determined  externally,  l.e.,  by  the  at- 
mosphere; this  coupling  is  clearly  most  effective  In  the  MAB  for  long 
periods  of  a week  or  more  (of.  Section  3). 

An  alternate  forcing  scheme  of  Interest  Is  that  of  an  Isolated  North 
Atlantic  storm  that  Impinges  on  the  MAB,  One  could  envision  an  Initially 
quiescent  SAB  and  then  consider  the  process  by  which  the  storm-induced 
MAB  vortlclty  disturbance  would  Influence  the  SAB.  From  the  discussion 
surrounding  Fig.  26,  all  of  the  significant  low  mode  energy  will  propagate 
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southwestward  In  the  MAB.  Upon  "rounding  the  corner"  at  Cape  Hatteras, 
however,  an  Increasingly  large  fraction  of  short  wave  energy  win  be  back- 
scattered  as  the  disturbance  proceeds  Into  the  SAB.  The  spectral  charac- 
teristics and  the  total  energy  of  the  disturbance  will  change  as  It  travels 
In  the  SAB,  with  long  wavelengths  (and  therefore  long  periods)  becoming 
relatively  more  important  than  short  wavelengths,  but  with  the  total  energy 
decreasing  due  to  the  short  wave  fraction  lost  to  backscatterlng.  The 
increasing  Importance  of  long  period  motions  as  one  proceeds  upstream  In 
the  SAB  is  suggested  in  the  interstation  longshelf  sea  level  coherence 
(Figs.  14  and  15)  and  In  the  MET-SL  coherence  (Figs.  16  and  21),  more 
notably  so  In  the  winter  case. 


5.  Summary  and  Conclusions 


The  principal  objective  of  this  study  was  to  see  if  continental  shelf 
waves  could  propagate  "around  the  corner"  at  Cape  Hatteras,  providing  a 
means  of  communication  between  the  Middle  and  South  Atlantic  Bights  (MAB 
and  SAB).  An  additional  objective  was  to  investigate  the  efficiency  of 
atmospheric  weather  disturbances  as  forcing  functions  for  these  waves.  The 
underlying  motivation  for  the  study  was  the  idea  that,  if  the  communication 
between  the  Bights  were  effective,  it  could  provide  a mechanism  for  intro- 
ducing vortlcity  disturbances  caused  by  North  Atlantic  storms  into  the 
Gulf  Stream.  The  data  sets  used  to  test  the  hypothesis  consisted  of 
coastal  tide  gage  and  meteorological  time  series  from  an  array  of  stations 
spanning  the  interbight  region. 

Analysis  of  the  tide  gage  data  showed  that  disturbances  with  periods 
longer  than  one  week  propagated  southwestward  in  both  Bights,  with  no 
apparent  interruption  at  Cape  Hatteras.  Similar  long  period  motions  have 
been  observed  as  far  north  as  Nantucket  Shoals  and  as  far  south  as  Key 
West,  so  it  is  probably  reasonable  to  consider  the  entire  eastern  United 
States  continental  margin  as  a single  "waveguide"  for  long  period  motions. 
The  wavelengths  implied  by  coastal  phase  lags  are  very  imprecisely  indicated 
by  this  study  to  be  of  the  order  of  1000  km,  which  is  approximately  the 
long  shelf  scale  of  the  observation  station  array  used  here. 

Several-day  period  motions  were  found  to  be  prominent  in  the  SAB  but 
evanescent  north  of  Cape  Hatteras  in  the  MAB.  Coastal  sea  level  phase  lags 
indicated  upstream  propagation  in  the  SAB,  but  the  lags  were  considerably 
smaller  than  would  be  expected  for  a monochromatic  several -day  period  shelf 
wave.  Also,  the  phase  lags  were  very  non-uniform  and  of  questionable 
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statistical  significance  over  most  of  the  SAB,  leading  to  the  conclusion 
that  a monochromatic  wave  is  an  Inappropriate  descriptor  of  the  several -day 
period  motions  In  the  SAB. 

Wind  stress  and  atmospheric  pressure  fluctuations  were  coupled  with 
the  sea  level  fluctuations  In  different  proportions,  depending  on  the  period 
of  the  motion,  the  season,  and  the  location  of  the  tide  gage.  Generally, 
both  wind  stress  components  (cross  shelf  and  longshelf)  provided  important 
forcing  for  periods  longer  than  about  a week  in  both  Bights,  and  atmospheric 
pressure  (after  correction  for  coherent  wind  stress  effects)  was  Important 
for  several -day  periods  in  the  SAB.  Strongest  coupling  occurred  In  the  SAB 
in  the  summer,  and  In  both  Bights  in  the  winter,  presumably  reflecting  sea- 
sonal differences  in  prevailing  weather  patterns. 

Consideration  of  the  energy  dispersion  characteristics  of  shelf  waves 
In  the  two  Bights  suggests  that  backscattering  of  long  period  motions  as 
they  enter  and  propagate  through  the  SAB  may  explain  the  presence  of 
several-day  period  energy  there.  This  result  is  a consequence  of  the  highly 
variable  bottom  slope  and  shelf  width  In  the  SAB,  and  Is  also  dependent  upon 
the  location  of  the  Gulf  Stream  relative  to  the  shelf  profile.  If  this 
hypothesis  is  correct,  the  absence  of  short-period .motions  In  the  MAB  Is  a 
result  of  the  relatively  uniform  shelf  profile  and  the  absence  of  the  Gulf 
Stream  there,  and  not  due  to  MAB-SAB  differences  In  the  atmospheric  forcing 
details.  It  Is  not  entirely  clear  that  scattering  Is  the  prime  source  of 
short  period  SAB  motions,  because  strong  and  very  selective  dynamical 

• 

coupling  between  atmospheric  pressure  and  sea  level  occurs  for  short  periods 
in  the  SAB.  The  latter  is  an  unexplained  result,  since  the  equations  Indi- 
cate that  direct  forcing  by  pressure  fluctuations  should  be  unimportant. 
Furthermore,  the  notion  that  pressure  correlation  with  sea  level  Is  just 
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an  indication  of  pressure  correlation  with  wind  stress  is  ruled  out,  at 
least  in  a linear  relationship,  because  the  pressure  transfer  functions 
upon1  which  the  above  statements  are  based  were  corrected  for  the 
pressure- to-wind  stress  linear  correlation.  A related  possibility  is  that 
the  pressure-to-sea  level  correlation  is  an  artifactual  consequence  of  the 
relationship  between  the  pressure  and  the  wind  stress  curl ; this  notion  is 
attractive  because  the  curl  is  more  likely  to  have  the  requisite  short 
period  energy  than  is  the  undifferentiated  wind  stress. 

It  is  clear  from  this  study  that  provincial  treatments  of  ocean  dynam- 
ics on  the  eastern  United  States  continental  margin  may  be  inappropriate, 
particularly  for  periods  of  a week  or  longer.  The  modeler's  field  of  view 
must  be  widened  sufficiently  to  encompass  the  different  but  coupled  dynamics 
of  the  two  Bights,  if  an  appropriate  understanding  of  the  low  frequency 
motions  is  to  be  gained.  The  dynamical  processes  that  seem  most  In  need  of 
treatment  in  light  of  these  results  are:  1)  scattering  and  energy  disper- 
sion of  shelf  waves  by  longshelf  coastal  and  topographic  irregularities; 

2)  statistical  details  of  the  process  by  which  realistic  wind  stress  and 
wind  stress  curl  fields  can  drive  oceanic  motions  over  the  shelf;  and  3) 
the  interaction  of  shelf  waves  with  the  Gulf  Stream,  particularly  with 
regard  to  momentum  and  energy  exchange  between  the  two  on  critical  surfaces 
where  the  wave  speed  vanishes  relative  to  the  mean  current. 


APPENDIX  A 

Time  series  of  low  pass  filtered  sea  level  and  meteorological  data.  See 
Section  2 for  details  of  filter  characteristics,  calendar  start  and  stop 
times,  and  wind  vector  coordinate  rotation  angles.  The  order  of  sea  level 
subset  presentation  follows  that  in  the  text,  proceeding  from  south  to 
north  along  the  coast. 
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2 Winter  adjusted  sea  level. 
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Fig.  A. 4 Winter  unadjusted  sea  level. 
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Fig.  A.6  Winter  unadjusted  sea  level. 
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Fig.  A. 8 Summer  adjusted  sea  level 


■ 

K 

EfKKmi 

IHifliMWJl 

nni 

kf  iwMAT1 1 

frBPmiiVi 

mmm 

mmji 

Elevation  (m) 


30  40  50  60  70  60  90  100  110  120 


Fig.  A. 


Time(days) 


10  Summer  unadjusted  sea  level. 
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Fig.  A.  12  Summer  unadjusted  sea  level. 
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Fig.  A. 15  VJinter  wind  vector  stick  plots.  The  ordinate  scale  gives 
the  vector  magnitude. 
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Fig.  A. 16  Summer  wind  vector  stick  plots.  The  ordinate  scale  gives 
the  vector  magnitude. 
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APPENDIX  B 


Multiple  coherence  and  partial  transfer  functions  from  the  regression  model. 
The  naming  convention  for  the  multiple  coherence  between  the  lettered  input 
meteorological  variables  and  the  output  sea  level  is:  a)  pressure  only, 
b)  pressure  and  cross  shelf  wind  stress,  c)  pressure,  cross  shelf,  and  long- 
shelf  wind  stress.  The  dashed  line  results  from  using  adjusted  sea  level 
as  the  output  for  case  c).  The  95%  null  hypothesis  level  for  case  c)  is 
about  0.45.  Transfer  function  phases  and  moduli  are  shown  for  both  wind 
stress  components  versus  unadjusted  sea  level  and  for  pressure  versus  unad- 
justed and  adjusted  sea  level.  The  error  bars  show  the  95%  confidence  inter- 
vals on  phase,  when  less  than  360  degrees.  The  regression  results  are  hier- 
archially  organized  by:  a)  season  (winter  then  summer),  b)  type  (coherence 
then  transfer  function),  and  c)  geographically  (north  to  south). 
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Fig.  Bl.  Winter  multiple  coherence. 
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Fig.  B3.  Winter  multiple  coherence 
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Fig.  B5.  Winter  multiple  coherence 
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Fig.  B8.  Winter  partial  transfer  functions. 


WILMINGTON  WIND  STRESS  VS.  GARDEN  CITY  SEALEVEL  CW> 


Fig.  BIO.  Winter  partial  transfer  functions. 
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Fig.  B17.  Summer  multiple  coherence 


9 


HATTER  AS  MET  t-OROLOGY  vs.  AVON  SEALEVEL  (S) 


adj.  sealvl 


FREQUENCY  ( C.  P.  D. ) 


PERIOD  (DAYS) 


Fig.  B18.  Summer  multiple  coherence 
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Fig.  B19.  Summer  multiple  coherence. 
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Fig.  B20.  Summer  multiple  coherence 
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APPENDIX  C 


Spectra  (m2/CPD),  coherence,  and  phase  (degrees)  for  sea  level  station  ad- 
jacent pairs,  proceeding  southward,  winter  and  sunnier  cases.  Negative 
phase  means  that  the  first-named  station  in  the  caption  box  leads  the 
second-named  station.  The  95%  null  hypothesis  level  is  shown  for  coherence 
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